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The leading twist transverse momentum dependent parton distribution function h 1T , which is 
sometimes called "pretzelosity," is studied. We review the theoretical properties of this function, and 
present bag model predictions. We observe an interesting relation valid in a large class of relativistic 
models: The difference between helicity and transversity distributions, which is often said to be a 
'measure of relativistic effects' in nucleon, is nothing but the pretzelosity distribution. Pretzelosity is 
chirally odd and can be accessed in combination with the Collins effect in semi-inclusive deep inelastic 
scattering, where it gives rise to an azimuthal single spin asymmetry proportional to sin(3<?!> — (f>s). 
We discuss the preliminary deuteron target data from COMPASS, on that observable and make 
predictions for future experiments on various targets at JLab, COMPASS and HERMES. 
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I. INTRODUCTION 



Processes like semi-inclusive deep-inelastic lepton nucleon scattering (SIDIS), hadron production in e + e~ annihi- 
lations or the Drell-Yan process [1-16] factorize at leading twist [17-19] allowing to access information on transverse 
momentum dependent fragmentation and parton distribution functions (TMDs) [20, 21]. The latter contain novel, 
so far unexplored information on the nucleon structure. In order to be sensitive to "intrinsic" transverse parton mo- 
menta it is necessary to measure adequate transverse momenta in the final state, for example, in SIDIS the transverse 
momenta of produced hadrons with respect to the virtual photon. 

The eight leading- twist TMDs /f, /f^?, gf, g^ T , hi, h±£, /ijy, h^ a [9], and further subleading- twist structures 
[22, 23] describe the structure of the nucleon in these reactions. The fragmentation of unpolarized hadrons is described 
in terms two fragmentation functions, D\ and H^ a , at leading- twist. In SIDIS (with polarized beams and/or targets, 
where necessary) it is possible to access information on the leading twist TMDs by measuring the angular distributions 
of produced hadrons. Some data on such processes are available [24-42]. 

The fragmentation functions and TMDs in SIDIS and other processes were subject to numerous studies in literature 
[43-68]. This is true especially for the prominent transversity distribution hi or the 'naively time-reversal-odd' 
functions like the Sivers function /f^?, the Boer-Mulders function hj~ a and the Collins fragmentation function H^ a . 
The so far probably least considered function is the 'pretzelosity' distribution hjji . 

The purpose of this note is therefore three-fold. First, we will review what is known about hj^. Second, we 
will calculate this function in the bag model. Third, we will present estimates for the transverse target single spin 
asymmetry proportional to sin(30 — 4>s) in SIDIS in which hj^S enters, and discuss the prospects to measure this 
asymmetry in experiments at Jefferson Lab. 



II. TMDS AND SIDIS 



Hard processes sensitive to parton transverse momenta like SIDIS are described in terms of light-front correlators 

dz~d 2 ZT 
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(1) 



We use light-cone coordinates a ± = (a ± a 3 )/^/2. In SIDIS the singled our 3-direction is along the momentum of 
the hard virtual photon, and transverse vectors like pr are perpendicular to it. The symbolically indicated gauge-link 
depends on the process [14-16]. In the nucleon rest frame the polarization vector S — (0, St, Sl) with S T + S\ = 1. 
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The information content of the correlator (1) is summarized by eight leading-twist TMDs [9], that can be projected 
out from the correlator (1) as follows 



1 

2 tr 



h 



E jk p> T S T 

M N 



At 



1 

2 tr 



7+75 <j>{x,pr) 



Sl9i 



Pt ■ St ± 



i 

- tr 
2 



9 J fc , o 4 h± {p T P T -\p^ k )S T 



T h 1 - 



e jk p T 
M N 



hi 



(2) 
(3) 
(4) 

= 1 and zero 

else. Integrating out transverse momenta in the correlator (1) leads to the three 'usual' parton distributions known 
from collinear kinematics jf(x) — J d 2 pT jf(x,p£) with j = f,g,h [85, 86]. Dirac-structures other than that in 
Eqs. (2, 3, 4) lead to subleading-twist terms [22, 23]. The fragmentation of unpolarized hadrons is described by a 
correlator analog to (1) parameterized in terms of two leading- twist fragmentation functions, D\ and Hi a . 

The process of SIDIS is sketched in Fig. 1. The cross section differential in the azimuthal angle <fih of the produced 
hadron, see Fig. 1, (and possibly differential also in other kinematic variables) has the following general decomposition 
[7, 69], where the dots indicate power suppressed terms, 



where the space-indices j, k refer to the plane transverse with respect to the light-cone and e 1 
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In F^y Sht the first index X — U(L) denotes the unpolarized beam (longitudinally polarized beam with helicity A). 
The second index Y = U(L,T) denotes the unpolarized target (longitudinally, transversely with respect to the virtual 
photon polarized target). The superscript reminds on the kind of angular distribution of the produced hadrons with 
no index indicating an isotrop ^-distribution, and ips is the azimuthal angle of the transversely polarized target. 

Each structure function arises from a different TMD. The chirally even /'s and g's enter the observables in connection 
with the unpolarized fragmentation function Df, the chirally odd /i's in connection with the chirally odd Collins 
fragmentation function Hi a . One has 
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HADRON PRODUCTION PLANE 




LEPTON SCATTERING PLANE 

FIG. 1: Kinematics of the SIDIS process IN — > I'hX and the definitions of azimuthal angles in the lab frame. 
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More precisely, TMDs and fragmentation functions enter the respective tree-level expressions in certain convolution 
integrals, indicated by <g> in (6-9), which entangle transverse parton momenta from TMDs and fragmentation functions. 
(Going beyond tree-level description requires introduction of soft factors [18, 19] from which we refrain here.) 

In general, such convolution integrals cannot be solved, unless one weights the DIS counts with an adequate power 
of transverse hadron momentum. For the sin(30/i — 0s) asymmetry that would mean to weight the events with P% ± . 
The resulting model-independent tree-level expression for that structure function is then given in terms of certain 
transverse moments of h^ T and [9] . 

The analysis of such single spin asymmetries (SSA), however, is involved because acceptance effects are difficult 
to control. So far only preliminary data on (some different) P^x-wcighted SSAs are available [32]. Moreover, the 
P/i_L-wcighting obscures the actual size of the effect. We shall therefore consider here the 'unweighted' SSA 



— sin(30-0s) 
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dependence with the known y-dependence removed, i.e. jr™( 3 ^ _ ^ s ) — (1 — y)F^y 3 s and Fuu = (l — y+y 2 /2)Fuu, 
such that the SSA depends only on x and z. (The SSA depends, of course, also on Q 2 . In the parton model this 
dependence is weak, logarithmic, and 'hidden' in the scale dependence of the distribution and fragmentation functions, 
which we do not indicate for brevity.) 

In this case, however, the integrals convoluting the transverse momenta can only be solved, if one assumes some 
model for the transverse momentum dependence. Here we shall assume the Gaussian Ansatz 

hj${x,ft) = hif{x) ex P(~Pf /P'v) > Ht a (z,K 2 ) = Ht a (z) e M-Kj/Kl) 

This is, of course, a crude approximation. However, besides being convenient [8], this Ansatz is also phenomenologically 
useful, provided the transverse hadron momenta are small compared to the relevant hard scale, (Ph±) <C Q in SIDIS, 
and one is interested in catching the gross features of the effects [52]. A high precision description of pp-effects requires 
methods along the QCD-based formalism of [2], see [70] and references therein for examples. 
With the Ansatz (11) we obtain for the SSA the following result 



d rin(3»-» fl ), „x _ C Gauss J2ge 2 a xh^ T 1)a (x)H^ {1/2)a (z) 
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which we have written such that the dependence on the Gauss model is 'dumped' into the factor 

3 

C G auss = f(z Pav /K av ) with /(a) = - a _ 2/3)3/2 • (13) 
In (12) we have introduced the transverse moments 

K& 1)a (x) = J d 2 p T Jjjj^ h&{x,tf) , Ht {1/2)a {z) = J d 2 K T iM Ht a {z,K 2 ) . (14) 

Little is known about the Gauss model parameters, p av and K w . Apart from constraints from positivity, we know 
nothing about p av and have only vague constraints on the Gaussian width of the Collins function [58, 59]. It is 
therefore gratifying to observe that /(a) in (13) is a slowly varying function of o satisfying < /(a) < 3/(2-^/2). This 
means that 

3 

< Ccauss < C" max = -^-j= . (15) 

For physically reasonable values of the ratio p av /K av (let us guess, for example, that it coincides within a factor of 
two with the corresponding ratio of the unpolarizcd widths p U np/^unp ~ 1.4 [55]) one may therefore expect that the 
factor Ccauss is within 20 % close to its maximum value. In the following, when we will be interested in estimating 
the maximal effects for the SSA (12), we will approximate Ccauss by C max in Eq. (15). 
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III. WHAT DO WE KNOW ABOUT hif ? 

Let us summarize briefly, what we already know, or will learn in the next sections, about the pretzelosity distribution. 

(1 ) It is chirally odd. 

(2) It has a probabilistic interpretation [8]. 

(3) At large- x it is predicted to be suppressed by (1 — x) 2 compared to /f (x) [65-67]. 

(4) It is expected to be suppressed also at small x compared to /f . 

(5) It must satisfy the positivity condition [68] 

< \ - < A>,pD (16) 

or, after integrating over transverse momenta: 

\h^ 1)a (x)\<l^(x)-g a 1 (x^j </?(x). (17) 

(6) Adding up (17) and the Soffcr inequality |/i?(x)| < |(/f + gf)(x) [71], we obtain the remarkable bound: 

\htW a (x)\ + \ht(x)\<ft(x). (18) 

(7) In the limit of a large number of colors N c in QCD pretzelosity is predicted, for pr ~ 0(N®) and xA c ~ 0(N®), 
to exhibit the following flavor dependence (the same relation holds for antiquarks) [72] 

\hf Wu (x) - ^ (1)d (x)| » |^ (1)u (x) + ^ (1)d (x)| (19) 
„ ' „ ' 

0(JV c 2 ) 0(N C ) 

(8) There are estimates from spectator model [73], and bag model, see Sec. IV. 

(9) In a large class of models, including the bag and spectator models, pretzelosity is the difference of helicity and 
transversity distributions, and in this sense a measure for relativistic effects, see Sees. IV and V. 

(10) In some sense it 'measures' the deviation of the 'nucleon shape' from a sphere [74]. 

(11) In simple (spectator- type) models, it has been related to chirally odd generalized parton distributions [75]. 

(12) It requires the presence of nucleon wave-function components with two units orbital momentum difference, e.g. 
s-d interference of quadratic in p-wave component [67] . 

(13) There is no gluon analog of pretzelosity. 

Some comments are in order. 

Concerning the point (4): This expectation is based on experience with evolution properties of integrated parton 
distributions. The evolution of the chiral-odd h\ differs significantly from that of the chiral-even /f or g\. If at some 
initial scale h\ were as large as /f or g\ (or the Soffer bound [71]) in the small- a; region, then at higher scales it would 
be suppressed [76, 77]. This pattern is expected also for other chirally odd TMDs. 

Concerning the point (6): For w-quarks, that usually dominate in the proton in the valence-x region, the bound 
(17) is more restrictive than the Soffer bound for transversity, because g\ (x) is positive. Thus, positivity gives more 
"room" for |/ii(x)| compared to \hf^ u (x)\ (and vice versa for the d-flavor). But from this observation alone we 

cannot draw any conclusion on the relative size of \h\(x)\ and \hf± 1)a (x)\, because neither function is a priori "forced" 
to saturate its allowed positivity bound. However, there are first indications that h\(x) is large, definitely larger than 

\{fi ~ 9i){x) [57-59]. So one may expect |^ 1)m (x)| < h\(x) on the basis of the first indications [57-59] . 

Concerning the point (7): The same flavor dependence as in Eq. (19) holds for all TMDs in a polarized nucleon. 
Notice, that there is also an unintegrated version of (19), see [72]. 

Concerning the point (13): Actually, in the decomposition of the gluon analog of the correlator (1) a structure 
appears that in Ref. [75] has been called hfj., which was of notational convenience for that work. (In Ref. [78] it was 
given a different name.) This gluon TMD, however, has different properties compared to our quark pretzelosity. For 
example, the of [75] is 'odd' under time-reversal while h^f with a = q, q is 'even'. 
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IV. PRETZELOSITY IN THE BAG MODEL 



In the MIT bag model, the quark field has the following general form [79-81], 

& a (x,t)= ^ N(nK){b a (nKm)ip nKjm (x,t) + dl(nKm)tp- n - Kjm (x,t)} , 

n>0,re=±l,m=±l/2 

where b' a and dJ a create quark and anti-quark excitations in the bag with the wave functions 

|x|- 
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For the lowest mode, we have n = 1, k = —1, and u>i t -i ~ 2.04 denoted as w = uj\ y -\ in the following. In the above 
equation, a is the 2x2 Pauli matrix, Xm the Pauli spinor, Rq the bag radius, x = x/\x\, and ji, are spherical Bessel 
functions. Taking the Fourier transformation, we have the momentum space wave function for the lowest mode, 



if m (k) = iV^KNRl 



*o(A:)xr, 



a-k ti{k)xr, 

where k = k/k with k = \k\ and the normalization factor TV is, 



(22) 
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(23) 



The two functions t i} i — 0, 1 are defined as 



U(k) — J u duji(ukRo)ji(ucj) . 
o 



(24) 



From the above equations, we see that the bag model wave function Eq. (22) contains both S and P wave components. 
Especially, to represents the S- wave component, whereas t\ represents the P-wave component of the proton wave 
functions. 

With the above wave functions, we can calculate all the leading order TMD quark distributions. Assuming SU (6) 
spin-flavor symmetry of the proton wave function, we have the up-quark and down-quark distributions, 



where h^ T can be written as, 



h^ T (x, kj_) = A 



fc 2 



A = 



P r , = - 



16a 



1 

3 ' 



7r 2 ( W -l)j 2 HM 2 



(25) 



(26) 



where ui is the lowest root of the bag eigen-equation as we mentioned above. The momenta k z and k are defined as 
k z = xMn — to/Ro, and k = \Jk 2 z + k\ . Ro and Mn are bag radius and proton mass, respectively. In our calculations, 
we fix the dimensionless parameter RqMn = Auj. Because h^ T depends on i 2 , clearly it is proportional to the squared 
of the P-wave [67], and is sensitive to the quark orbital angular momentum in the proton. 

In order to study the positivity bounds we have discussed in the last section, we list the unpolarized and polarized 
quark distributions in the MIT bag model. We have 



fi(x,k±) 
gi(x,k±) 
hi(x,k±) 
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where according to the SU(6) spin flavor symmetry 

fl(x,k ± ) = N q f 1 (x,k±), N u = 2, N d = l, (30) 

and the polarized TMDs of definite flavor, g\ and h\, are given by the expressions analog to (25). 
From the above results, we read off the relations, for example 

fi(x,k±)+g 1 {x,k±)= 2h 1 (x,k ± ) (31) 

hi{x,k 1 _)-h^ 1 \x,k 1 _) = fi{x,k A _) (32) 

where hfj^\x,kj_) = k\/{2M] <[ )hi T {x,k±). Thus, in the bag model out of the four functions fi, gi, hi, h^ T only 
two are linearly independent. 

In QCD the various TMDs are all independent of each other, and describe different aspects of the nucleon structure. 
In simple models, however, it is natural to encounter relations among TMDs. Since they are all expressed in terms 
of to and ti representing the S and P-wave components of the proton wave function, Eqs. (22, 24), the TMDs cannot 
be all linearly independent of each other. The specific form of the relations (31, 32) can be traced back to Melosh 
rotations, which relate longitudinal and transverse nucleon polarization states in a Lorentz-invariant way [82]. The 
relations (31, 32) arise therefore in some sense from Lorentz invariance properly took into account in a simplified 
description of the nucleon. 

Considering the flavor factors (25, 30) dictated by SU(6) symmetry, we obtain from (31, 32) the relations 

f!(x, k±) + g\{x, k ± ) = 2h\{x, k ± ) (33) 

h\{x, k ± ) - h^ 1)q (x, k ± ) = ^ fi(x, k±) . (34) 

The latter can be rewritten as follows (notice that the signs are such that the first part really remains an equality) 

\hl(x,k ± )\ + \h^ 1)q (x,k ± )\ = \Mf*( x , k± ) < ft(x,k±) . (35) 

Comparing (35) with (18) we see that the inequality (18) is always satisfied in the bag model, and actually never 
becomes an equality. In fact, transversity and pretzelosity explore only | (only |) of their allowed positivity bound 
for the it-flavor (rf-flavor). Also the inequality (16) is satisfied, because the Soffer inequality holds in the bag model, 
which follows analogously from (33), cf. also the detailed discussion in [83]. 

Concerning large- N c , we recall that in the quark model formulated for a general (odd) number of colors N c , the 
flavor factors (25, 30) are given by P u = (N c + 5) /6 and P d = (-N c + 1)/6 while N u = (N c + 1)/2 and N d = (N c - 1)/2 
[84]. Thus, the bag model results respect large- N c predictions (19) because the employed SU(6) symmetry does so. 
We also see, however, that l/iV c -corrections are sizeable for N c = 3. 

Since only two out of the 4 functions fi, gi, hi, h^ T are linearly independent there are numerous relations among 
them. The probably most interesting relation results from subtracting the relations (31, 32) from each other. This 
yields (since here a common flavor factor, P q , enters) the result 

g\{x, k±) - h\{x, k±) = hjfi )q (x, k±) . (36) 

This is result is remarkable from the point of view of the popular statement 'the difference between transversity and 
helicity distributions is a measure for relativistic effects in nucleon' [86]. Here, in the framework of the bag model, we 
can quantify this statement: The difference of helicity and transversity is (a transverse moment of) pretzelosity. 

From this discussion we immediately see what happens to h^jl in the non-relativistic limit: it vanishes. In this 
sense, we arrive at the conclusion that the pretzelosity distribution itself is a measure for relativistic effects in nucleon. 

Let us turn to the discussion of the numerical results 1 . Fig. 2 show results for the quantity 



h^(x) — J d 2 kj^h^(x,kj^) 



(37) 



The MIT bag gives also rise to antiquark distributions, though to unphysical ones as the ff(x) come out negative, violating positivity. 
Nevertheless these contributions are crucial for the normalizations f dx ff(x) = N q or the momentum sum rule J2 q I dxxff(x) = 1, 
that are satisfied only if one includes the contribution from negative x (where ff(x) means minus the antiquark-distribution). Moreover, 
the TMDs receive non- vanishing support also from \x\ > 1, and the sum rules are satisfied only when integrating over the whole rr-axis. 
It has been discussed in literature how to deal with these caveats, see for example [87] . In this work, we limit ourselves to quark TMDs 
at < x < 1 (we do not discuss 'valence distributions' /j val (a;) = fii x ) — fi( x ))- When discussing sum rules, however, we shall integrate 
the TMDs over the whole rr-axis. 



7 



Bag model results for scale-dependent quantities refer to a low initial scale Qq which is understood to be of the order 
of a typical hadronic scale, say Q = 0(M N ), or even lower, see e.g. [88]. 

The u-quark distribution is negative. Compared to the d-quark distribution, it has opposite sign to and its absolute 
value is four times larger — according to (25). The pretzelosity distribution functions have opposite signs compared to 
the transversity distribution functions. Interestingly, the h v ^{x) are larger than the transversity distribution functions 
h\(x) in the bag model, and also larger than ff{x). However, it has to be noted that this quantity, as defined in 
Eq. (37), is not constrained by positivity bounds. The h-^ix) reach minima or maxima, depending on the flavor, 
around x = (0.2-0.3). 

What is constrained by positivity are the transverse moments hf^ 9 (x) of the pretzelosity distribution, see 
Eqs. (16, 17, 18). Our results for these quantities, which are defined in Eq. (14), are shown in Fig. 3. The transverse 
moments of the pretzelosity distributions are smaller than the transversity distributions. For x > 0.1 their absolute 
values are always more than a factor 3 smaller than those of h\ (x) . The extrema of the transverse moments are shifted 
towards larger x compared to h^(x), and appear around x — (0.3-0.4). 



hjj(x) 




0.2 0.4 0.6 0.8 x 



h 1T (x) 




0.2 0.4 0.6 0.8 



FIG. 2: The parton distribution function hf^,(x) vs. x from the bag model (results obtained here) in comparison to f1(x) and 
h\{x) from the same model. The functions hf^(x) are rather large, even larger than ff(x). Notice, however, that h^(x) itself, 
as defined in (37), is not constrained by positivity bounds. All results refer to the low scale of the bag model. 
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FIG. 3: The transverse moment hf^ q (x) vs. x from the bag model (results obtained here) in comparison to fi(x) and h\(x) 
from the same model (rescaled by the factor ±| for better visibility). The functions h^f (x) satisfy the positivity bounds, 
Eqs. (17, 18). The results refer to the low scale of the bag model. 
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V. PRETZELOSITY IN THE SPECTATOR MODEL OF [73] 



The pretzelosity distribution was calculated in the spectator model [73] . Strictly speaking valence quark distribu- 
tions were computed there. However, in that model calculation the difference between 'quarks' and 'valence quarks' 
can be disregarded to a good approximation. In the spectator model of Ref. [73] pretzelosity is given by 



hi${x, k±) = J2 Lr ' q h ^{x, k±) , = | , xj l / u = i ; ifl* = o , L a ' d = 1 



(38) 



where the sum goes over contributions from a scalar spectator diquark (s), and an axial vector spectator diquark (a). 
The other TMDs have the same 'decompositions' (38), which is dictated by the SU(6) symmetry. The contributions 
of the various diquark spectators to pretzelosity, and for comparison, the TMDs /i, g\ and h\ read [73] 



fi(x,k±) = J r (a;,fc_L) 

gl(x,k±) = a r J r (x,k±) 

hi(x,k±) = a r J r (x,k±) 

hif(x,kj_) = a r J r (x,k±) 



(xM N + to) 2 + k]_ 
(xM N + to) 2 - k\ 



(xM N + to) 2 



-2M 2 N 



where the 'spin factor' a r assumes the values a s = 1 and a a = —\, and 
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2(27T)3(fc 2 + ' 



A 2 = A 2 (l - x) + xMf. - x{\ - x)M% . 



(39) 
(40) 
(41) 
(42) 

(43) 



The meaning and the values of the parameters are: spectator masses M a — 0.8 GeV and M s — 0.6 GeV, quark mass 
to = 0.36 GeV, cutoff A = 0.5 GeV, and a = 2. The normalization factor N is such that J dx f{{x) = 1, and for the 
nucleon mass the rounded value Mm = 0.94 GeV was used [73]. Other choices of parameters are also possible, see [73] 
for details. The results refer to a scale as low as the initial scale of the GRV parameterization /Xq = 0.23 GeV 2 [89], 
and possibly even lower [73] . 

The contributions from the respective diquarks satisfy the relations, for example, 



a r f{(x, kj_) + g[(x, kj_) = 2h\{x, kj_) 
h[(x, k±_) - hf^ r (x, k±) = a r f{{x, k±) 



(44) 
(45) 



which arc the analogs of the bag model relations (31, 32) discussed in Sec. IV. But, since M a ^ M s , the relations 
(44, 45) in general do not give rise to relations among TMDs for specific flavors. More precisely, one recovers the 
relation (34) in the spectator model for d-flavor, since here the scalar quark decouples. For the it-flavor the relation 
(34) follows only in the limit M a — > M s , because in that limit the decompositions (38) coincide with the bag model 
SU(6) flavor relations (25, 30). Therefore, although both are based on the SU(6) symmetry, the spectator model is 
more general than the bag model with respect to the flavor dependence. 

Notice that in spectator models the mass difference between scalar and axial-vector diquarks is responsible for the 
A-nucleon mass splitting whose physical value is reproduced for M a — M s = 200 MeV [90] . In the large- N c limit we 
have Ma — Mjy = 0(N l T 1 ), which implies that also M a — M s — > in that limit. As for M a — M s the flavor relations 
(25, 30), are recovered, see above, the same discussion applies here as in Sec. IV. Thus, the spectator model respects 
the large- N c relation (19) — formally and practically, although with non-negligible l/V c -corrections, see below. 

From this discussion we learn an interesting lesson about the model dependence of relations among different (in 
QCD independent!) TMDs. Relations that involve unpolarized and polarized TMDs are 'flavor dependent', see 
(33, 34), and valid only in 'flavor-blind' models like the bag model, see Sec. IV and [83], or constituent quark models, 
see [91]. Such relations break down in general in models with non-trivial flavor structure. In fact, the different masses 
of scalar and axial-vector diquarks in the spectator model spoil such relations. 

What about relations involving polarized TMDs only? If we subtract the relations (44, 45) from each other f[{x, k±) 
drops out, and since in that model the flavor structure is precisely the same for all polarized TMDs (independently of 
the scalar and axial- vector diquark masses), we recover Eq. (36). Thus, the relation (36) among hclicity, transversity 
and pretzelosity seems to be valid in a larger class of models. We shall come back to this point in Sec. VI. 
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Let us now discuss numerical results. Fig. 4 shows h 1 ^,{x) defined in (37). The u-quark distribution is negative and 
(in magnitude) much larger than the positive d-quark distribution. The h^{x) are larger than h\(x) in that model, 
and also larger than f1(x) (we recall that this quantity is not constrained by positivity). 

Fig. 5a shows the transverse moments hf-f (x). One can convince oneself that the positivity constraints (17, 18) 
are always satisfied in the model of [73] independently of the specific values for the spectator masses. However, here 
the flavor structure is more involved, such that it is not possible to see a priori what fraction of the positivity bound 
the transverse moment of pretzelosity explores — in contrast to bag model, Eq. (35) — and it is worth to look at this 
in more detail. Fig. 5a shows that the transverse moments of pretzelosity never exceed even half of the (17). Also 
the positivity condition (18) is satisfied, see Fig. 5b. Noteworthy, pretzelosity and transversity use, especially for the 
it-flavor, a large fraction of the room allowed by positivity conditions. 

Next we confront the results from the bag model (obtained here, see Sec. IV) and the spectator model of Ref. [73]. 
There is a good qualitative agreement, the signs and magnitudes of hf^(x) and h^' q {x) agree, see Fig. 6a and b. 
The transverse moments h^S (x) are much smaller than the transversity distributions h\(x). With the exception of 




0.2 0.4 0.6 0.8 x 0.2 0.4 0.6 0.8 x 



FIG. 4: The parton distribution function h 1 ^,{x) vs. x from the spectator model of Ref. [73] in comparison to ff{x) and h\ (x) 
from the same model. The functions h^(x) are rather large, even larger than f1(x). Notice, however, that h 1 j,{x) itself, as 
defined in (37), is not constrained by positivity bounds. All results refer to the low scale of the model [73]. 



hjj ^ q (x) vs. bounds (a) 
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(b) 




0.2 0.4 0.6 0.8 x 



FIG. 5: (a) The transverse moment h^f (x) vs. x (solid lines). It does not exceed even half of the bound <J (a;) = i(/f ~ g\)(x) 

(dashed lines) allowed by positivity, see Eq. (17). (b) The combination \h^f, (x)\ + |/i*(x)| as function of x (solid lines). This 
combination must not exceed fi(x) (dashed lines), see Eq. (18), which is the case. 
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very small- a; they are at least 3 times smaller than h\ (x), see Fig. 6c. 

Let us finally compare the flavor combinations (hf^ u ± hfj, )(x) in Fig. 6d. As we have discussed, above and 
in Sec. IV, both models are conceptually in agreement with large- N c . But in both models, of course, the finite value 
N c = 3 is used. It is therefore of interest to investigate to which extent the models reflect the large- N c pattern in 
Eq. (19). In fact, Fig. 6d shows that the modulus of the (u — d) flavor combination is always larger than that of the 
(it + d) flavor combination, as one would expect on the basis of large- N c arguments in Eq. (19). But the difference 
between the magnitudes of the (it ± d) flavor combinations is not pronounced. This is true especially in the bag 
model where the 'large' (it — d) flavor combination is only |-times larger than the 'small' (u + d) flavor combination. 
(Recall that this factor is JVc 3 +2 , i.e. in fact large for N c — > oo.) The spectator model reflects more clearly the large- N c 
behavior predicted in Eq. (19) — especially if we recall that these predictions are valid for x ~ 0(1/N C ) [72]. 

To which extent can one expect models formulated for finite N c — 3 to respect large- iV c predictions of the kind 
(19)? Fig. 6d gives a flavor on that. In nature, large-iV c results are found to hold within a similar accuracy [92]. 
Interestingly, the Sivers function predicted to exhibit in the large- N c limit a flavor dependence analog to (19), see 
[72], seems to obey large- N c predictions rather closely [53-57, 60, 61]. 




hj ( T )q (x)/h q j(x) ( c ) hjy q (x) vs. large N c ^ 




0.2 0.4 0.6 0.8 x 0.2 0.4 0.6 0.8 x 



FIG. 6: Comparison of results from bag model (computed here, cf. Sec. IV) and spectator model, Ref. [73]. (a) h 1 ^,(x) vs. x. 
(b) hf± 1)q (x) VS. X. (c) The ratio of h^ 1)q (x) to hl(x) vs. x. (d) The flavor combinations (hf^ 1 ± h^ 1 ' ld )(x) vs. x. All results 
refer to the low scales of these models. 
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VI. HOW GENERAL IS THE RELATION IN EQ. (36) 

In the bag and spectator model the relation (36) is valid not only for quark but also for antiquark distributions. 
Hereby one has to take into account that the Dirac-structure in the chirally odd hi and h^ T is C-odd, while it is 

C-even for g\. Thus, for antiquarks the relation (36) reads g\ {x, k±) + h"(x, k±) = —hf^ a (x, k±) with a = u, d. 

While strictly speaking in the bag model anti-quark distributions arc unphysical, see footnote 1, one obtains in 
principle consistent results in the spectator model. In [73] no antiquark distributions were considered, however, this 
can be done precisely as sketched in Eqs. (38-43) but with additional parameters for the (four-quark-) spectators, 
which must be tuned. The results, however, are analog to (38-43) and satisfy the relation (36) for antiquarks. 

So, in the quark models we have discussed the relation (36) holds for both, quark and antiquark TMDs. It will also 
be valid in a large class of relativistic models, e.g. the constituent quark models of the kind used in [91], or in the 
model of [93] (which remains to be verified by direct calculations). The question that naturally arises is: How general 
is the relation (36)? 

When discussing this question it is important to observe that Eq. (36) relates chirally even and chirally odd TMDs, 
i.e. it relates quark (and anti-quark) TMDs having gluon 'partners' and such having no gluon counterparts, cf. Sec. III. 
It is clear that the evolution properties of these TMDs are different. Therefore, if the relation (36) were valid at some 
initial scale, then it certainly would not be valid at a different scale. 

These considerations imply, that the relation (36) can be valid only in 'no-gluon models'. This expectation is 
supported by the model calculations of Ref. [75]. There (among others) quark and gluon TMDs in a hypothetical 
'quark target' were computed. It was found that gi(x, k±) — hi(x, k±) ^ (meaning that the model is 'relativistic') 
but difference was not related to pretzelosity. The latter is zero in that model (to the considered order of 'a s ') [75]. 
This is in line with our expectations: the explicit inclusion of gluon degrees of freedom spoils (36). 

Nevertheless, Eq. (36) could turn out to be approximately satisfied, and useful for estimating pretzelosity on the 
basis of the helicity and transversity distributions, which we know presently better, or are on a good way to that 
[57-59]. At this point, since this relation is found in models without gluons (where we simply could 'neglect' the 
gauge-links), presumably the fcj_-integrated version of (36) is a more reliable prediction of relativistic quark models. 



VII. PRELIMINARY COMPASS DATA, AND FUTURE EXPERIMENTS AT JLAB 

In the COMPASS experiment the sin(30 — 0s) and other SSAs were measured on a deuteron target [41]. We shall 
estimate the deuteron distribution functions by neglecting nuclear binding effects and exploring isospin symmetry as 
(analog for u and d, and we neglect strange and heavier quarks): 



,-Lu/D _ , -Lli/p ,-Lu/n 

1 T — 1 T I < h i ^y 1 
,-Ld/D _ h ±d/p h ±d/n 



+ hiS ■ (46) 



Let us denote by N(ir + ,D) the numerator of the SSA A s ^ 3 ^ ^ in production of positive pions from the deuteron 
target. It is given by (we skip prefactors irrelevant for the qualitative discussion) 

N(ir+, D) = (hif + hi4){AHi im + H^ uni ) + (h^ + h^){H^ v + 4H^ uni ) (47) 

We estimate the maximum effect for this SSA as follows (and analogously for other pions) 

\N(tt+, D)\ < \h$ + h$\ ■ \4H^ av + Hi unf \ + \h$ + hi£\ ■ |#i Lfav + 4tf 1 - Lunf | (48) 

where (and analogously for antiquarks) 

\hf^ u + h^ 1)d \ < \hf^ u \ + \h^ 1)d \ < + ft g» - g() . (49) 

Using information on Collins effect [58, 59], the parameterizations [89, 94] for fi(x), gi(x), Di(z) at a Q 2 — 2.5 GeV 2 , 
and assuming that positive (negative) hadrons at COMPASS are mainly positive (negative) pions (this is legitimate 
to a good approximation), one obtains the results shown in Fig. 7. We compare the positivity-bound result to the 
preliminary data [41]. 

At larger x > 0.1 the statistics analyzed in [41] is not sufficient to be sensitive to the pretzelosity distribution. 
However, below x < 0.1 COMPASS data favor that pretzelosity does not saturate its bound. (This would be in line 
with our models. But those refer to low scales, and are strictly speaking not applicable to small-x.) 
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^™( 3< M>s)( x ) on deuteron (a) A^x^ '^(x) on deuteron (b) 
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FIG. 7: The transverse target SSA A*^ 3 ^ for deuteron estimated on the basis of the positivity bound vs. preliminary 
COMPASS data [41]. 

This could hint at the following. Either, the chiral-odd pretzelosity is suppressed at small- a; with respect to the 
chiral-even fx, which is expected, see Sec. III. Or, the u and d-flavors have opposite signs, as predicted in the large- N c 
limit and the models. Or, both. In view of the early stage of art it is too early to draw definite conclusions. 

An important observation, however, is that in spite of the small (consistent with zero) effect observed at COMPASS, 
pretzelosity does not need to be small in the region x > 0.1, see Fig. 7, where JLab can measure with great precision. 
So there is no discouragement due to small COMPASS pretzelosity effect for the planned CLAS measurements! 

It remains to be mentioned that HERMES also studied the sin(3</> — <ps) azimuthal modulation in the transverse 
target experiments. It was quoted that this SSA is zero within error bars, see e.g. [29, 31, 40]. Presumably, an SSA 
as large as what we obtain from saturating positivity (see below) would have been seen at HERMES. However, since 
no HERMES data were shown on that observable, we cannot draw quantitative conclusions from that. 

In order to see how further analyses from HERMES, as well as future experiments at COMPASS and Jefferson Lab 
could improve our understanding of pretzelosity, let us estimate the upper bounds for the modulus of the SSA for 
charged pion production from various targets. Using the positivity bound (17) we obtain the results in Fig. 8. 

The SSAs for charged pions could reach up to ~ 5%. Even asymmetries reaching half or one third of that size 
would be measured, especially in the region of x ~ (0.1-0.4) in the CLAS experiment. This is shown in Fig. 9, where 
we plot the ir + SSA from a proton target in the kinematics of CLAS with 12 GeV beam upgrade. Shown are also 
error projections for 2000 hours run time from [96]. Notice that the models predict a negative SSA. 

We add that the ir° SSA is small on any target, because the modulus of this SSA is proportional to \Hi v + \ 
which is small (actually zero within error bars), since ps — jFff Lfav holds within error bars [57-59]. 




FIG. 8: The modulus of the transverse target SSA /lyy 3 * < * s - ) in charged pion production from proton, deuteron, and neutron 
targets as function of x. Estimates on the basis of the positivity bound. 
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FIG. 9: The SSA A 



sin(3^>-0 s ) 
UT 



in 7r + production from a proton target in the kinematics of CLAS 12 as function of x. The error 
projections are from [96]. The shaded areas indicate the range allowed by positivity bounds. The bag model (obtained here) 
and spectator model [73] predict a negative sign for the SSA, i.e. in the lower shaded area. 



VIII. CONCLUSIONS 



We reviewed and discussed the properties of the pretzelosity distribution function /ij~ T , and presented a study of 
this leading- twist, chiral-odd, transverse parton momentum dependent distribution function in the bag model, and 
supplemented our findings with a detailed comparison to spectator model results [73] . 

In these models we observed an interesting relation, which we expect to be valid at low scales in a wide class of 
relativistic models. It can be summarized for illustrative purposes by the following assertion: 



helicity — transversity = pretzelosity. 



(50) 



That the difference between the helicity and transversity distributions is 'a measure of relativistic effects' is known 
since long ago [86] (and was also recognized in a bag model calculation). However, now we are in a position to 
make this statement more precise. This difference is just pretzelosity. Thus, h^ T 'measures' relativistic effects in the 
nucleon, and vanishes in the non-relativistic limit where helicity and transversity distributions become equal. 

This relation is not supported in models with explicit gluon degrees of freedom [75], and, of course, cannot be true 
in QCD where all (eight) transverse momentum dependent parton distribution functions are linearly independent. 
Nevertheless, the relation (50), see Eq. (36) for its precise formulation, could turn out to be a useful approximation. 
In view of the numerous novel functions involved, any well- motivated approximation is welcome and valuable [64] . 

Besides being useful for extending our intuition on relativistic spin-orbit effects in nucleon [67, 74], the relation (50) 
has also an important consequence on transversity. In the bag and spectator model hj^ is negative. Since g\{x) is 
positive, this implies that hf(x) > gi(x). For the d- flavor signs are reversed, but in any case |ft.f(a;)| > |fl?(aO|- This 
is found also in models, e.g. [95]. 

In the bag model, the negative sign of arises because it is proportional to minus the square of the p-wave com- 
ponent of the nucleon wave function. Thus, in models with no higher orbital momentum (d-wave, etc.) components, 
hjjt is manifestly negative {h^jt has opposite sign dictated by SU(6) symmetry, and predicted in large iV c [72]). 

This prediction can be tested at JLab. Since the production of positive pions from a proton target is dominated by 
the u-fiavor, one expects a negative sin(3</> — <ps) SSA, see Fig. 9. 

Forthcoming analyses and experiments at COMPASS, HERMES and JLab [96-98] will provide valuable information 
on the pretzelosity distribution function, and deepen our understanding of the nucleon structure. 
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